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ABSTRACT 



Properties are developed for the failure rate function 
and residual mixing probabilities resulting from a mixture 
of distributions. After a detailed investigation of these 
properties for purely exponential mixtures, two more general 
mixtures are considered. Emphasis is placed on relating the 
current failure rate in the mixed population to the popula- 
tion's residual composition. The results have potential 
application in determining an appropriate period of "burn in. 
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I. INTRODUCTION 



Suppose that at time t=0 there exists a population of 
similar elements. This population could be made up of a type 
of transistor, a type of hydraulic check valve, a certain type 
of lock, e;^ cetera . But whatever the elements composing the 
population may be, they are all designed to perform the same 
task. And being a functional physical object, each element 
will have some random lifetime, depending on a known or unknown 
probability distribution, during which it functions properly. 
When this random lifetime has elapsed, the element fails. 

The usual definition of the probabilistic failure rate 
function of an element with random lifetime will be used 
throughout this thesis. Namely, the failure rate function of 
an element is the density function for .the lifetime of the 
element divided by the element's survival function. The 
scenario can nov; be further developed by requiring that the 
population consist of subgroups of elements, v;ith each sub- 
group having its own unique failure rate function applicable 
to each element in the subgroup. There is, then, a mixture 
of failure rates in the population. It will be assumed that 
at time t = 0 the exact proportion of the total population 
held by each subgroup is known either explicitly in the dls-r 
Crete case or implicitly through a cumulative distribution 
function in the continuous case. And so, the proportion of 
the population at time t = 0 having a particular failure rate 
is knov/n. 
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Now that the population has been defined. It Is possible 
to ask the two questions whose answers this thesis Is con- 
cerned with. First, If an element Is selected at random from 
the population at time t = 0, what Is the nature of the 
randomly selected element’s failure rate function? And given 
that a particular element selected at random from the popula- 
tion at time t = 0 Is still functioning at time t > 0, what • 

Is the probability that that element was a member of a 
specific subgroup? 

For someone wanting an explanation of the practical Impor- 
tance of these questions, a plausible example might be the 
following. Suppose that the population consists of transistors 
which were not marked as to subgroup membership due to some 
oversight during the production process. If the Installation 
of one of the transistors Into a piece of electronic equipment 
at time t = 0 Is required, what can be predicted about the 
failure rate of the Installed transistor? And on the basis 
of Its performance, what can be said about that transistor's 
subgroup source? 

Others have already developed some of the answer to the 
first question raised above about the nature of the failure 
rate of an element selected at random from the population at 
time t = 0. For Instance, Barlow, Marshall, and Proschan 
[Ref. 1] have shown that a mixture of distributions each 
having a nonincreasing failure rate Itself has a nonincreasing 
failure rate. Proschan [Ref. 2] applies this result to help 
explain an observed phenomenon In the aircraft Industry and 
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to conclude that a mixture of exponentials will have a 
decreasing failure rate. Esary, Marshall, and Proschan 
[Ref. 3] show that a mixture of distributions with decreasing 
hazard rate average has a decreasing hazard rate average. 

And Gnedenko e^ a^. [Ref. 4] indicate that in a mixture of 
exponentials the failure rate always approaches as a limit, 
V7lth Increasing time, the least exponential parameter value 
of all elements composing the population. This thesis sheds 
further light on the nature of the failure rate of an element 
drawn at random from a population composed of elements with a 
mixture of failure rates. No claim is made that the exposi- 
tion here begins to give a complete answer, and areas needing 
further study are indicated. 

The second question, about the conditional probability 
that a randomly selected element has come from a subgroup 
with a particular failure rate function given that the element 
is still functioning at a certain time, does not seem to have 
been given much attention in the literature. And once again, 
the results presented here are certainly not exhaustive. More 
research is required if the question is to be completely 
answered. 

Briefly, the procedure follov/ed in this thesis is to show 
results for increasingly more complex mixtures of exponentials. 
Once this has been done, some attention is given to more 
general mixtures. 
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II . THE CASE OF TV/0 EXPONENTIALS (DISCRETE PARAMETERS) 



Let the population at time t = 0 consist of a proportion 
Pj^ of elements distributed exponentially with parameter 
and a proportion P 2 of elements distributed exponentially 
with parameter The equality p^ + P 2 = 1 holds. Assume 

that is less than X 2 * It is a well known result that the 
failure rate functions h^^(t) and h 2 (t) for elements in the 
two respective population subgroups are given by h^(t) = 
and h 2 (t) = X 2 (see Figure 1). 

At time t = 0 an element is selected at random from the 
population. The following definitions are made: 

Aj^ is the event that the element selected comes from 
the subgroup with failure rate X^. is the event 

that the element selected comes from the subgroup 
with failure rate X 2 * T is the lifetime of the 
element selected at random from the population. And 
F^(t), f^(t), and h^(t) are the element’s survival, 
density, and failure rate functions, respectively. 
Using these definitions, some equations can be derived: 

P^(t) = P(T>t) = P(T > t lA^)P(A^) + P(T > t 1A2)P(A2) 



= 



+ P2e"^2t. 



f^(t) = -F^’(t) = p^X^e + p2X2e"^2t ^ 
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Figure 1 
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(Xi - X2)t 



hrp(t) = 






p^x^^ + P2X26 



+ ?2 ® 



(Xi - X2)t 



hrp' (t ) = 



-p^P2(X^ - X2>^ ■ ^2)t 






+ p,e 



(Xi - X2)t' 



h,p’*(t) = p^P2(X2 - X^)-^e 



3^(Xi - X2)t 



Pt - PoS 



(Xi - X2)t 



[Pl 



+ P2© 



(Xi - X2)t 



r. 



Using the definition of h^(t), the following theorem 
results . 

Theorem 1. For 0 < p^ < 1, the relationship < h^(0) < X2 holds ♦ 
Proof. Note that 

h,p(0) = p^X^ + P 2 X 2 = (1 - P2^^i ?2^2 ^1 P2^^2 “ ^1^* ^ 

The following theorem is a special case of the general 
result shown in Ref, 1. 

Theorem 2. The failure rate h^(t) ^ a decreasing function . 
Proof. Prom the expression for h^'(t) given above it is seen 
that h^' (t) <0 for every t. # 

Reference indicates the general result of which the 
following is a special case. 

Theorem 3- lim h,j^(t) = X^. 
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Proof. The result follows directly from the expression for 
h^(t) and the assumption made above that < X 2 » ^ 

Theorem 3 is an interesting result since it means that, 
no matter how great the probability P 2 of choosing an element 
at time t = 0 from the subgroup with failure rate X 2 may have 
been, the failure rate h^(t) eventually becomes as close to 
X^ as desired. An intuitive explanation of this phenomenon 
is that as the element chosen at random from the population 
is observed to have an ever increasing lifetime without failure 
it becomes ever more likely that it is an element with a low 
failure rate. But the lowest failure rate it can have is X^ 
since there were only two possible failure rates in the popu- 
lation from which the element was taken, and X^ was the least 
of these two. Intuitively, it would seem that the longer the 
element functions the more likely should its failure rate be 
X^. The limiting process of Theorem 3 reflects this likelihood. 

The range of h^(t) is bounded. 

Theorem 4. X^ <h^(t) < X 2 for every t. 

Proof. The proof follows from combining the results of 
Theorems 1, 2, and 3* ^ 

Some general statements can be made about the shape of the 
failure rate curve. 

Theorem 5 . For p^ < P 2 j the function h^(t ) has an inflection 
point . Before the time at v;hlch this inflection point occurs 
the function is concave, and it is convex after that time . 
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Proof. The function h^(t) has an inflection point when its 
second derivative vanishes. 

- p,e(^ - = 0«t = in g) . 

For <P2> this is a feasible time. That is, it is not a 
negative time. 

Since p-j^<P2> it is seen from the expression for h,j,'’(t) 
that this function has a negative value at t = 0, and it is 
monotone increasing for t >0. Therefore, h,j(t) is concave 
before the time at which its inflection point occurs, and it 
is convex afterwards. § 

The value of the failure rate function at the time the 

Xl + 

inflection point occurs is given by h_(t^) = , where 

- X^) (p2y ' 

As a result of Theorem 5 the graph of h,j,(t) for the general 
case of a mixture of two exponential failure rates has the 
form shovm in Figure 2. An intuitive explanation for the 
shape of this curve can be given. Prior to the time at which 
the inflection point occurs h,j,(t) always has values closer to 
X2 than x^. It seems that if the element chosen at random 
has a relatively brief lifetime it probably has a high 
failure rate. The concavity of the curve, its shape tendency 
towards X2, prior to the inflection point seems to be an 
Indicator that the element is more likely to have come from 
the subgroup with failure rate X2. The failure rate X^ 
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Failure Rate Function Resulting From a 
Mixture of Two Exponential Types (p^ < P2) 

Figure 2 
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seems to have a stronger attraction for the curve than 
does prior to the inflection point. 

At the inflection point, the failure rate function h,p(t) 
has been shown to have a value exactly midway between and 
X^- The intuitive explanation for the inflection point 
occurring when this happens is that the element is probably 
equally likely to have come from either subgroup. The 
failure rates X^ and X2 seem to have an equal attraction for 
the curve. 

After the time the inflection point occurs, the curve 
becomes convex. The intuitive explanation for this is 
similar to that given above for the concave portion of the 
curve . 

For the case when ?2 Pi» ^ip’ ' ('t) is zero only at a 
negative time. And so for all feasible times, t^O, the 
function h^(t) is convex and closer to X^ in value than to 
X2* The intuitive explanation for this Involves the same 
reasoning used above for the case p^ < P2* The case p^ = 

P2 = ^ is the same except that h^(t) has an inflection point 
at t = 0 , and it then becomes convex. 

A precise understanding about the conditional probability 
that an element randomly selected has come from a particular 
one of the two subgroups given that the element is still 
functioning at a certain time requires the definition of a 
residual mixing probability, called a "mixing probability" 
hereafter. Though the definition in this section is for the 
tv/o exponentials case, it will in later sections have an 
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easy extension to more complex mixtures. By the mixing 
probability at time t is meant the function Pj^(t), i = 1 or 2, 
such that Pj.(t) = P(A^|T>t), t^O. This may be rewritten as 



P^(t) = P(Aj^|T > t) 



P(A., T>t) 
P(T > t ) 



P(T > t|A^)P(A^) 
P(T > t ) 



Pj_e 



-Xn- 1 






For the particular case of 1 = 1, this means that 



P^(t) = 



e-^1^ + P2e"'2t 



Pi + 



P2^ 



(Xi - X2)t 



Similarly , 



2 P2 Pi® ^ ^ 



Note that P^CO) = p^ and P2(0) = P 2 * 

Once the mixing probabilities P^(t) and P 2 (t) have been 
expressed in the above forms, it is easy to see that the 
intuitive explanation given above for the shape of the h^(t) 
curve is correct. Since X^ < X 2 , Pj^(t) is monotone increasing 
to the value 1, and P 2 (t) is monotone decreasing to a value 
of 0. And for the case where h^(t) has an inflection point 
(the case when p-j^<P 2 )j it can be seen that the following 
relationships hold: 
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Pj^Ct) < P2(t) If t < tj 



P^(t) = P2(t) if t = tj 
P^(t) > P2(t) if t > tj . 

The mixing probabilities can be used in another way to 
see how the curve h,p(t) arises. 

Pl^l^ + P2^2® 

h^(t) = = X^P^(t) + A2P2(t) = E^(A). 

-Xit , -Xot 
p^e J- + P20 ^ 



Here A is the random variable taking values X^ and X 2 at 
time t with probabilities P^(t) and P 2 (t), respectively. So 
once these conditional probabilities are known, the function 
hip(t) is determined as an expected value. This is another 
fact having extension to cases more complex than a mixture 
of tv7o exponentials, as will be seen in later sections. 

The first derivative h^'(t) can be represented in terms 
of expected values at any time t too. 



dP^(t) 

dt 



2 2 
( E p.e"^<3^)(-X.p.e-^i^) + (p I X.p.e"^jh 

.1=1 ^ iJ: I j=l ^ ^ 

( E p.e"^J^)^ 

J=1 



= {E,(A) - X,}P.(t) for 1 = 1, 2. 
0 11 
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And so 



h^’(t) = E^’(A) = + X2P2’(t) = -{E^(A^) - E^^(A)} 

= - Var^ ( A ) . 

This is another way of showing that h^p'Ct) < 0. 

The second derivative h^’’(t) also can be put in a 
concise expected value form. 



h,p'’(t) = -{E^»(a2) - 2E^(A)E^' (A)} 

2 2 2 
= -{ I X/P '(t) - 2E.(A) E X.P '(t)} 
i=l 1 ^ i = l ^ 

2 2 

= - Z - 2X.E. (A)}P. • (t) 

1=1 ^ 1 r 1 

2 2' 

_ E {X. - 2X.E. (A)HX. - E.(A)}P.(t) 

** ^ 1 u 1 w 1 

= E^{[A^ - 2AE^(A)][A - E^(A)]} 

= E^{A^ - 3A^E^(A) + 2AE^^(A)} 

= E^{A - E^(A)}^ . 



The expression of h,p’’(t) as the third central moment of A 
will also have extension to general mixtures of exponentials. 
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III. MIXTURE OF THREE EXPONENTIAL TYPES (DISCRETE PARAMETERS) 



Before considering in detail the general case of a mix- 
ture of n exponential types, some results which have been 
found for a population composed of three exponential types 
are presented. Function definitions for the n exponentials 
case are given first, and then in the remainder of this 
section it is assumed that n = 3- The more general defini- 
tions are needed in the next section. The definition of A. , 
i = 1, 2, ..., n, parallels that given in the last section 
for and Also it will be assumed that the failure 

rate functions of the population subgroups satisfy 



_ n 

P„(t) = P(T>t) = E P(T > t|A.)P(A. ) 

^ i=l 1 ^ 



? --Xft 
E p . e 1 

1=1 ^ 



f^Ct) = 



_ n 

- Frn' (t ) = E p. X. e 
i=l ^ 



-X j[t 



E p . X . e 1 
i=l ^ ^ 



hi(t) ^ 



E p . e 
1=1 ^ 



-Xj[t 



When the definition of the mixing probabilities P^(t) 
given in the last section is extended to populations where 
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the number of subgroups is n>2, the failure rate function 
can once again be written as an expected value. 



n 

h^^Ct) = Z X.P. (t) = E. (A), 
i i=i 1 1 ^ 



Using this fact and changing the upper index on all summation 
signs in the last section from 2 to n, it is apparent that in 
the general case 



h,p' (t ) = -Var^(A) . 

h^' ' (t) = E^{A - E^(A)}^. 



A program was written in FORTRAN and run on the IBM 360 
as a means of observing the behavior of h,p(t) during various 
time intervals. For n = 3j examples checked with the aid of 
the program led fo a conjecture which has been proven to be 
true . 

Theorem 6. ^ the population is composed of a mixture of 

three exponentials ^ with ^2 ~ ^1 ^3 ~ ^1 ^ 

and p^ then h^(t) has an inflection point at 

tj = ^ In^^j . And h^(t^) = X^. 

Proof. Note that at time t^, v/hich is feasible since Pj^_^Po, 



PlCtj) = P^(tj). 



Therefore, 

h^(tj) = = (^2 " a)P^(tj) + a)P^(tj) 
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And also 



h^”(t) = E^{A - E^(A)}^ 



= {(X 2 - a) - X 2 >^P^(t) + (X 2 -X 2 )^P 2 (t) 



+ { (X 2 + a) - ^ 



= - P^(t)} = 0. 
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Theorem 6 details the case for n = 3 which is analogous 
to the n = 2 case, where it was shown in the last section 
that the inflection point, if there is one, occurs at that 
time t when h^(t) is exactly midway between X^ and X 2 * 
Theorem 6 describes the case for n = 3 when the inflection 
point will occur exactly midway between X^ and X^. And at 
the time this inflection point occurs, for both the n = 3 
and n = 2 cases, the least and greatest Indexed mixing 
probabilities are equal. Even the forms for the two cases 
of the times at which this occurs are analogous. For n = 2 
it was shown that, for P-]^£P 2 j ^he inflection point occurs 
when 




Theorem 6 says the inflection point occurs when 
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It would be nice to find that h^(t) always has an inflec- 
tion point when it takes on a value midway between and X^. 
But this result does not carry over to more complex mixtures, 
and it falls to hold even for n = 3 when not all of the 
conditions of Theorem 6 are satisfied. As an example, 
consider the case where 

X^ — 1 p^ — . 2 

X^ = 2 P 2 = .2 

X3 = *1 P3 = .6 . 

Here 2.5 is midway between X^ and X^, but the Inflection 
point was found to occur when .24 < t < .25 and 2.5739 < h^(t) 

< 2.5915. 

Theorem 6 does not extend to the case n = 4, for instance, 
even vjhen conditions analogous to those required in the 
theorem are satisfied. As an example, 

X^— 1 p^— .1 

X 2 = 2 P 2 = .2 

^3 = 3 p^ = .3 

X|^ - 4 P2| = . 4 . 

The Inflection point was found to occur when .4469 < t < .4470 

and 2.50955 < h^(t) < 2 . 50960 . 

A particular property applicable here of a general result 
proven in the next section is that P^(t) approaches the value 
1 while P 2 (t) and P 3 (t) approach 0 as t becomes infinite. So 
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just as for the case n = 2, the least indexed mixing 
probability has 1 as a limit with increasing time, and all 
other mixing probabilities have 0 as a limit. 
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IV. GENERAL MIXTURE OF EXPONENTIAL TYPES (DISCRETE PARAMETERS) 



The general case in which the original population con- 
sists of n subgroups, with all elements in any subgroup 
having the same exponential failure rate differing from the 
failure rate associated with every other subgroup, has some 
interesting properties. Definitions of p^ and P^(t) carry 
over directly to the general case from the n = 2 and n = 3 
cases. In the last section expressions were given for P^(t), 
hfp(t), h;p'(t), and h,p’’(t) for this mixture of n exponential 
types. 

Theorem 3 is just a particular case of a property pos- 
sessed by a general mixture of exponentials. 

Theorem 7. 11m hfp(t) = X^. 



Proof. 



h^(t) 



n 

P "|^1 ^ 

1=2 ^ 



(Xi - Xl)t 



n 

p., + I p.e 
1=2 ^ 



(Xi - Xi)t 






Reference 1 shows that the function h,p(t) is monotone 
decreasing. This of course means that h,p’(t) <0 for all t. 

When n = 2, h^(t) has at most one inflection point. 

This result has not yet been proven for the general case, 
even though computer runs of examples with various values of 
n have all revealed at most one inflection point for h,p(t). 
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If it could be proven that h^**(t) is monotone Increasing, 
then it would follow immediately that h^(t) has at most one 
inflection point. But this does not hold in general. 
Consider, for instance, an example where n = 5 and the param- 
eter values for the mixture are 



Here h^**(t) is negative at t = 0, increases to a positive 
maximum at t = .824, and then begins to decrease. All that 
can be concluded is that if h^(t) is concave at t = 0 it 
must have at least one inflection point because h^(t) 
approaches as a limit as t becomes infinite. 

A study of mixing probabilities in the general case 
reveals a number of interesting properties. Writing P^(t) 
in the form 




Pj^(t) 




it can be seen that 



lim P . (t ) 
t-»-® ^ 
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This is an intuitively plausible result in light of the 
limiting behavior of h^(t) discussed above. The longer the 
element selected at time t = 0 continues to function, the 
closer its failure rate h,p(t) gets to It is reasonable, 

then, to expect that the conditional probability that the 
element selected came from the subgroup with failure rate 
X^, given continuing survival, should become greater with 
increasing time. And since the h^(t) curve, for large t 
at least, continually moves farther from the failure rates 
^2’ •••» reasonable that the conditional proba- 

bilities that the element selected has come from a subgroup 
with a failure rate greater than X^ should become less and 
less. 

Additional facts about the behavior of the mixing proba- 
bilities can be determined. 

Theorem 8. ^ h^(0) > X^, Pj^(t) reaches a maximum when 

h,j(t) = X^. 

Proof. An Immediate extension to the general case of a 
result shown for n = 2 is 

dP, (t) 

= {h^(t) - x^)P^(t). 

This means that 
dP. (t) 

— ^ > 0 h^(t) > X^ 

dP. (t) 

— d t ' — = ° ^ 
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i 



dP^(t) 



0 ^ h^(t) < . 



dt 



< 






The last line of the above proof shows that If h,p(0) < X^ 
the function Pj^(t) has its maximum at t = 0 and is monotone 
decreasing thereafter since h,p(t) is monotone decreasing. 

The intuitive explanation for Theorem 8 seems to follow 
from the monotone property of h^(t). Since h,p(t) is always 
decreasing, it must, if h,p(0) > X^, get closer and closer to 
assuming the value X^ until the time when h^(t) = X^. It 
can be expected, then, that the conditional probability the 
element chosen came from the subgroup with failure rate X^, 
given survival until time t, should be increasing to a maxi- 
mum at the point in time when h,p(t) = X^. As time continues 
beyond this, hrj,(t) gets farther and farther from X^. It 
would seem it must be getting less likely that the element 
chosen came from the subgroup whose elements have failure 
rate X^^. 

A property expressed in the proof of Theorem 6 for the 
case n = 3 does not carry over to an arbitrary exponential 
mix. In Theorem 6 it was shown that there is a case for 
n = 3 such that when h^(t) takes on a value midway between 
X^ and X^ the total mixing probability mass for X^ above 
this midvalue is equal to the total mixing probability mass 
belov; the midvalue. That is to say, when h^(t) reaches the 
midvalue, the element is equally likely to have come from a 
subgroup with failure rate greater than the midvalue as it 
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Is to have come from one with failure rate less than the 
midvalue. But the example for n = 5 used earlier in this 
section shows this result does not carry over to the general 
case, even though the X^'s form an arithmetic progression 
just as the X^'s in Theorem 6 did. In the example, h,j(t) = 

X^ at t = .38276. But mixing probabilities at this time for 
X^ values below the midvalue X^ satisfy P^(t) + P 2 (t) = .42433, 
while mixing probabilities for X^ values above the midvalue 
satisfy P^jCt) + P^(t) = .45834. 
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V. CONTINUOUS PARAMETER EXPONENTIAL MIXTURES 



So far, the only populations which have been considered 
are those v;ith a discrete number of subgroups and, thus, a 
discrete number of parameters. In this section it will 
again be assumed that each element in the population has 
exponential lifetime indexed by a failure rate parameter, 
but now the parameter will not be restricted to discrete 
values only. Instead, it will be assumed that it is a con- 
tinuous, nonnegative random variable A with cumulative 
distribution function Gy^(x) and density g^(x). When 
is a step function, with G.(X ) = 1 for finite X , the result 
is one of the cases already considered. 

Definitions in this continuous case are extensions of 
the ones in discrete cases. For instance. 



F^(t) = P(T>t) = / P(T>tlA = X)g^(x)dx = / e"^^g (X)dx, 

a a 

where a > 0 and Gy^(a) = 0. 

A number of the general properties of h^(t) are the same 

here as vjhen the parameters were discrete. Reference 1 shows 

that hrp(t) is monotone decreasing. Reference 4 indicates 

that lim h™(t) is the minimum value that the random variable 
t-»-» ^ 

A can assume. And at time t = 0, h^(t) is the expected 
value of the random variable A. 
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Theorem 9* h,p(0) = E(A). 

Proof. Note that 



F^(t) 

h,p(t ) 



oo 

/ e"^^dG^(X) = E(e"^^) . 



-F^’(t) = E(Ae . 



E(Ae~^^) 

E(e"^^) 



And therefore 



h^(0) = E(A) 



§ 



The properties of hrp(t) are illustrated in the following 
example. 

Example. Suppose that A is uniformly distributed over the 
interval [a, b], a > 0. Then 

^A^^^ ^ b - a * 



b -xt,, 
a 



-at -bt 
e - e 

(b - a)t 



rrj.it) 



(1 -I- ab)e~^^ - (1 + bt)e 
(b - •a)t^ 



w _ (1 + at)e"^^ - (1 + bt)e 

T^^^ 7"^::at -bt. , 

(e - e )t 
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Using elementary limiting arguments, it is found that 



h„(0) = 11m h„(t) = 

^ t->0 ^ 




And also 

11m hm(t) = a . 

t-»-oo 

The analog in the continuous case of mixing probabilities 
is a conditional density function, g^(X). This function 
gives the conditional density of the random variable A, 
given that the element selected at time t = 0 is still 
functioning at time t. For and such that 

AX = X. , - X. is small, {g. (X)}AX will, then, give a good 
approximation to the probability that an element surviving 
until time t is one with failure rate A satisfying 
X^ £ A ^ ^1+1' conditional density is defined by 



Note that, as in discrete cases, the failure rate func- 
tion can be interpreted as an expected value at any time t 
since 





00 



e“^^g^(u)du 



/ e 
a 



f xe ^^dG^(X) 



A 



00 





= f Xg^(x)dX = E(AlT>t) = E^(A). 
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Theorem 8, too, has its continuous analog. 

Theorem 10. ^ h^(0)>X, g^(X) reaches a maximum when 

h^(t) = X. 

Proof. Note that 



-ut 



-Xt 



-Xt 



-ut 



dg (X) ® g(u)du][-Xe g(X)] + [e g(X)][/ ue g(u)du] 



dt 



00 2 
[/ e“'^^g(u)du] 



= {h^(t) - X}g^(X) . § 

And of course, as before, if h,p(0)_<x the function g^(^) 
has its maximum at t = 0 and is monotone decreasing there- 
after since h^(t) is monotone decreasing. 
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VI. SOME NONEXPONENTIAL MIXTURES 



In this section the requirement that all elements in 
the population have exponential failure rate is relaxed. 
Although a completely general mixture is not considered, 
previous theory is extended to mixtures more complex than 
the exponential ones considered thus far. 

A useful function in studying more complex mixtures is 
the hazard function. This function has been used by others 
(see, for instance. Ref. 3)- The hazard function of an 
element with random lifetime is the function H(t) defined 
by H(t) = -In F(t), where P(t) is, as usual, the element's 
survival function. And if f(t) and h(t) are the element's 
density and failure rate functions, respectively, the follow- 
ing relationships hold. 

h(t) = = H' (t) . 

P(t) 

f(t) = -P'(t) = H'(t)e“^^^^ . 

And since 

t 

H(t) = / h(s)ds , 

0 

it follov7s that 

f(t) = h(t)e"^^^^ . 

Note that H(t) is an increasing function since P(t) is 
decreasing. 
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Now suppose that the population is composed of n subgroups. 
All elements in subgroup i have density function f^(t), 
survival function P^(t), failure rate function h^(t), and 
hazard function H^(t). None of these functions for a sub- 
group is equal to its corresponding function in any other 
subgroup. And again at time t = 0 each subgroup i makes up 
a proportion p^ of the total population. Without any other 
restrictions, this would describe the most general mixture 
of a discrete number of subgroups. 

Consider now a case where each failure rate h^(t), i > 1, 
is greater than h^(t) by a constant at all times t. That is, 

h^(t) = h^(t) + 

h2(t) = h^(t) +k2 

hn(t) = h^(t) + k^ ,■ 

where 0 = k^ < k^ < ... < Then 

t 

H. (t) = / h.(s)ds = H, (t) + k.t . 

0 ^ 1 X 

And since h^(t) < h^Ct) < ... < h^(t), H^(t) < H 2 (t)< ...<H^(t). 
And also 

f^(t) = h^(t)e"^i^^^ = [h^(t) + . 

P. (t) = . 
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An element selected at random from the population at 
time t = 0 will, then, have random lifetime described by the 
failure rate function 



z p . f . (t ) 

Z p.F (t) 
i=l ^ 



n 

Z 

1=1 

n 

Z 

i=l 



p^[h^(t) + 

p.e"[^Hi(t) + kit] 



And since F^(0) = 0 by assumption, it follows that 
n 

h„(0) = Z p.h (0) . 

^ 1=1 ^ ^ 

This corresponds to the result shown for exponential cases 
that hip(0) is the expected value, weighted by the p^, of 
failure rates in the population at time t = 0. In fact, just 
as before, h^(t) can be interpreted at all times t as an 
expected value since now the mixing probabilities are given 
by 






p.F. (t) 
^11 



I P.F.(t) 

j=l ^ ^ 



p g-[H^(t) + kit] 

I ' 

3=1 •> 



And so 



h,p(t ) 



n 

z 

i=l 



h^(t)P^(t) . 



3 ^ 



Pjj^(t) can actually be put into the same form that the 
mixing probability had in the exponential case since 



P^(t) 



? -kit 

Z p . e J 

j=l ^ 



This means that the limiting result shown for exponential 
cases holds here. That is, 

'"l for i = 1 

lim P. (t) = 

t-»-<» ^ [o for i = 2, n. 

Another limiting result holds too. If h,p(t) is put in 
the form 



h^(t) = h^(t) + 



n 

Z 

i=2 



p.kie-'^i* 



n 

Z 

i=l 



p^e 



-ki t 



5 



it can be seen that lim [h„(t) - h,(t)] = 0. Just as in the 

t-j-oo t 

exponential case, h,j,(t) ultimately approaches as a limit the 
least failure rate function present in the original population . 

Theorem 8, too, has its counterpart. 

Theorem 11. ^ h^(0) > h^(0), P^(t) reaches a maximum when 

h^(t) = h^(t). Otherwise , P^^Ct) ^ continually decreasing . 



35 



Proof. The argument exactly parallels that used before, 
since 



dP^(t) 

dt 



{h^(t) - hj^(t)}P^(t) . 



And h|p(t) gets continually closer to h^(t) since 



Z e 

i=2 ^ ^ 



n 

z 

i=l 




t 



is a decreasing function. 



# 



Consider now a population in which the failure rates for 
subgroups satisfy 

h^(t) = k^h^(t) 
h^Ct) = k2h^(t) 

hn(t) = kj^h^(t) , 

v/here 1 = k^ < k 2 < ... < Note that in this case 

F^(t) = [F^(t)]^i. 

The result in the exponential case that h^(t) is a 
decreasing function cannot be expected to carry over at all 
times here since h^(t) may at some times be Increasing. 
However, h^(t) can be expressed as a time dependent multiple 
y(t) of h^(t), and y(t) is a decreasing function. To see 
this it is helpful to consider the case where the multipliers 
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are not discrete k.'s but instead come from a continuous 

1 

random variable K ^ 1 with cumulative distribution function 

Gj^(k). When Gj^(k) is a step function, and Gj^(k^) = 1 for 

some finite k , the discrete case is described, 
n’ 



F^(t) = / [P^(t)]^dGj^(k) . 

00 

f,r(t) = / k[F^(t)]^”^f^(t)dGj^(k) . 

h^(t) = h^(t)y(t) , where the function y(t) is 

00 

/ k[F^(t)] dGj^(k) e[KF, (t)^] 

y(t) = i_ = K - • 

/ [F (t)]^dG„(k) E[F^(t) ] 

1 1 K 



Note that 

y'(t) {/ [F^(t)]^dGj^(k) } ^ 

00 oo 

= / [F^(t)]^dGj^(k) / k^[F^(t)]^"^[-f3^(t)]dGj^(k) 

oo 00 

- / k[F^(t)]^dGj^(k) / k[F^(t)]^"^[-f3_(t)]dGj^(k). 

This means that y*(t) < 0 if and only if 

^ / k[F^(t)]^dGj^(k)}^ < / [F^(t )]^dGj^(k) / k^[F^ (t ) ]^dGj^(k) . 
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by the Cauchy-Schwarz inequality for Stieltjes integrals. 

This proves that y(t) is a decreasing function. And this 
means that once the hfp(t) curve crosses below an h^^^Ct) curve, 
it never intersects that hj,(t) curve again. This result, 
too, was found in the exponentials mixture case. 

Returning to the discrete mixture, h^(t) again approaches 
the least failure rate function present in the population 
as time Increases. 

Theorem 12. 



Proof. Using the definition of hrp(t), it can be seen that 






1=2 



k 



1 



1 . 



+ r - •'l)Hl(t) 



i=2 



# 



Noting that 



it can be seen that h,p(t) has the usual expected value 
Interpretation at any time t. That is, 

n 

h^(t) = h^(t)P^(t) = E^[hj^(t)] . 

Theorem 11 holds here also since once again 

dP. (t) 

■ ■ - = (h^(t) - h. (t)}P^(t) . 

More research is needed to determine the existence of 
possible necessary and sufficient conditions which failure 
rate functions in the original population need to satisfy 
in order that properties of h^(t) and Pj^(t) seen for mixtures 
considered in this thesis hold. 
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